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Abstract 

We investigate the heavy quark mass effects in the virtual photon structure func- 
tions F^{x, Q 2 , P 2 ) and F2{x, Q 2 , P 2 ) in the framework of the mass-independent renor- 
malization group equation (RGE). We study a formalism in which the heavy quark mass 
effects are treated based on parton picture as well as on the operator product expan- 
sion (OPE), and perform the numerical evaluation of F^ s (x, Q 2 , P 2 ) to the next-leading 
order (NLO) in QCD. 
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§1. Introduction 



The Large Hadron Collider (LHC) 1 ' has started its operation and it is anticipated that 
the signals for the new physics beyond the Standard Model (SM) will be discovered. Once 
these signals are observed, more precise measurements will need to be carried out at the 
future e + e~ collider, so-called the International Linear Collider (ILC)P In such cases, it 
is still important for us to have detailed knowledge of the SM predictions at high energies 
based on QCD. 

It is well known that, in high energy e + e~ collision experiments, the cross section of 
the two-photon processes e + e~ — > e + e~ + hadrons dominates over that of the one-photon 
annihilation processes e + e~ — > 7* — > hadrons. The two-photon processes provide a good 
testing ground for studying the predictions of QCD at high energies. Here we consider the 
two-photon processes in the double-tag events where both of the outgoing e + and e~ are 
detected (see Fig. [T|). In particular, we investigate the kinematical region in which one of the 
photons with momentum q is far off-shell (large Q 2 = — q 2 ) while the other with momentum 
p is close to the mass-shell (small P 2 = — p 2 ), can be viewed as a deep- inelastic scattering 
where the target is a photon rather than a nucleonP In this deep-inelastic scattering off 
photon targets, we can study the photon structure functions, which are the analogs of the 
nucleon structure functions. 




e 



e 



Fig. 1. Deep-inelastic scattering on a virtual photon in the e + e collider experiments 

The unpolarized (spin-averaged) photon structure functions F^ix, Q 2 ) as well as F^(x, Q 2 ) 
of the real photon (P 2 = 0) were first studied in the parton-model (PM)® and then investi- 
gated in perturbative QCD (pQCD). The leading order (LO) QCD contributions to F^ix, Q 2 ) 
and F2(x, Q 2 ) were obtained by WittenP The next-to-leading order (NLO) corrections to 
F^ix, Q 2 ) were calculated by Bardeen and BurasP These results were obtained in the 
framework based on the operator product expansion (OPE)0 and the renormalization group 
equation (RGE)P The same results were rederived by the QCD improved PM.^» 
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The structure function F%{x, Q 2 , P 2 ) for the case of a virtual photon target (P 2 ^ 0) was 
investigated in the LO and in the NLO by pQCDp^^ Also the virtual photon structure 
function F2(x,Q 2 ,P 2 ) was studied in the LOP- 1 ' In fact, these structure functions were 
analyzed in the kinematical region, A 2 < P 2 < Q 2 , where A is the QCD scale parameter. 
The advantage of studying a virtual photon target in this kinematical region is that we 
can calculate the whole structure function, its shape and magnitude, by the perturbative 
method. This is contrasted with the case of the real photon target where in the NLO there 
exist nonperturbative pieces. Recently the QCD analysis was made for F^x, Q 2 , P 2 ) up 
to the next-to-next-to-leading order (NNLO) and for F2(x,Q 2 ,P 2 ) up to the NLO. 121 And 
more recently the target effects on these structure functions were studied^ and compared 
with the existing experimental dataP^'^ In these calculations all the relevant quarks were 
assumed to be massless. 

In this paper we examine the heavy quark mass effects on the photon structure functions 
F2(x,Q 2 ,P 2 ) and F2{x,Q 2 , P 2 ). Indeed, the heavy quark mass effects for the two-photon 
processes, especially in the deep-inelastic kinematical region, have been studied by many 
authors ^3$ MbBSiMB g u ^- they were not treated within the framework of the OPE and the 
RGE. Our analysis here is performed in the framework of the QCD improved PM powered 
by the parton evolution equations and based on the mass-independent renormalization group 
approach in which the RGE parameters, i.e., (3 and 7 functions, are the same as those of the 
massless quark case. We consider the system which consists of n/ — 1 massless quarks and 
one heavy quark together with gluons and photons. Then, the heavy quark mass effects are 
included in the RGE inputs; the coefficient functions and the operator matrix elements. In 
the case of the nucleon target, the heavy quark mass effects were studied by a method based 
on the OPE in Ref. I2T1) . where the heavy quark was treated such that it was radiatively 
generated and absent in the intrinsic quark components of the nucleon. This picture does 
not hold for the case of virtual photon target, since the heavy quark is also generated from 
the virtual photon target with light quarks at high energies. We should consider both the 
heavy and light quarks equally as the partonic components inside the virtual photon. 

In the next section, we derive the evolution equations for the parton distribution func- 
tions in the case where rif — 1 light quarks and one heavy quark are present. In section 3, we 
calculate the heavy quark effects in the virtual photon structure functions, F^(x, Q 2 , P 2 ), 
F2(x, Q 2 , P 2 ) and F^ s (x, Q 2 , P 2 ), and compare our theoretical predictions with the existing 
experimental data. The final section is devoted to the conclusions. We discuss the diagonal- 
ization of the anomalous dimension in Appendix A, and the parton-model derivation of the 
master formula in Appendix B. 
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§2. The evolution equations in the presence of the heavy quark effects 

In this section we consider the evolution equations for the case in which rif — 1 massless 
quarks and one heavy quark exist. The extension to the system with many heavy quarks is 
straightforward, since we can repeat this treatment recursively. 

2.1. Heavy quark effects and operator mixing 

We discuss the heavy quark effects in the RGE mixing and derive the master formula 
for the moments with quark mass effects in this subsection. Although the OPE is the useful 
formalism, we can get more physically intuitive picture in the PM approach, so we discuss 
the heavy quark effects in the evolution equations in the parton model. Now we consider 
the evolution equations for n/ — 1 massless quark parton distribution functions (PDFs) 
q l (x,Q 2 ,P 2 ) (i = 1, • • • ,rif — 1) and one heavy quark PDF q H (x,Q 2 ,P 2 ) together with 
the gluon PDF Q 2 , P 2 ) and photon PDF P 7 (x, Q 2 , P 2 ), Experimentally, this situation 

corresponds to system of u,d,s (massless quarks) + c (heavy quark) for kinematical region of 
PLUTO datsPJ and u,d,s,c (massless quarks) + b (heavy) for that of L3p^ We write down 
the DGLAP equations for q\ q H , G 1 , P 7 ; 



dq^Q^P 2 ^ 
rflnQ 2 



x dy 

i y 



+P qG [-, Q 2 j GT(y, Q 2 , P 2 ) + P i7 I -, Q 2 ) H(y, Q 2 , P 
dq H {x ) Q 2 ) P 2 ) 



g Pa (j, Q 2 ) g j (y, Q\ P 2 ) + Ph (~ g 2 ) g H (y, Q\ P 2 



ding 2 

*dy 

i y 



E p *> (j . Q 2 ) ^ Q 2 > p2 ) + p ** (j. Q 2 ) « H (y> Q 2 > p2 



+Phg (jj, Q 2 j CT{v, Q 2 , P 2 ) + Ph, (j, Q 2 ) I ' ( !)■ Q 2 . P 2 ' 

dG\x, Q 2 , P 2 ) 
dlnQ 2 

^E p G* (~> Q 2 ) i j (y> ^ p2 ) + p oE (X g 2 ) i H (y, Q 2 , p 2 ) 



X 



dy 

l y 



-Pgg ( g 2 J ^(y, g 2 , p 2 ) + p G7 f i g 2 ) / - 1 ,,. a 2 , p- 



X 



(2-1) 



[2-2) 



;2-3) 
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where .Pjj = 8ijP qq + — P£ is the splitting functions of j-parton into i-parton, and the first 



term represents the process that j-quark splits into i-quark without through gluon, and the 
second term stands for the splitting through gluon and P qq and P® are both independent of 



quark flavors, % and j. P qq is relevant for the flavor-singlet part and starts in the order of a 2 . 



We now define the singlet combination, q^, and the non-singlet part of i-th quark q l NS as 
well as the non-singlet combination q]^ s for the light-flavors as, 



E 

i=i 



Ins 



71/ - 1 



7 



9ns 



(24) 



i=i 



Note that Y^i=\ Ins = 0- We should remember that the photon PDF in the virtual photon 
P 7 does not evolve within the order c^qed, so we set P 7 (y, Q 2 , P 2 ) = 5(1 — y). We can rewrite 
the Eqs. fl2-ip . (I2 2|) and (12-31) in terms of q1,q]j S - This can be done by the following steps. 
At first, we sum up Eqs. (12- ip . (12-21) and (12-31) from i — 1 to i — Uf — 1 and we get the equation 
for the parton distribution defined as a row vector q 7 (x, Q 2 , P 2 ) = (ql, q 11 , G 7 , q N s), 



d 



,q\^Q 2 ,P 2 



dy 

y 



dlnQ 2 

where the splitting functions are given by a matrix 



P 



q\y ) Q 2 ) P 2 )P[- 1 Q 2 

y 



X 



+ k(x,Q 2 ,P 2 



( P s 

qq 


Plh 


Plg 
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Phl 


Phh 


Phg 





Pgl 


Pgh 


Pgg 





\ 








pNS 

m 1 



with 



p s 



HL 



P„ + 
-Ip 



n 



f 



1 ~. 



P* P, 



n 



s 

qq, 



Phh 



qq: 



LH 



1 ~ 



P q % P LG = (n f -l)P c 



qG, 



Pqq + —Pf 



qq- 



P 



HG 



P 



HG, 



P 



GL 



pS 



P 



GH 



P 



GH, 



P 



GG 



P 



GG, 



P 



NS 

qq 



qq 



(2-5) 



(2-6) 



(2-7) 



and the inhomogeneous term, k = (ki, ku, k G , kNs) describing the parton-photon mixing is 
given by 



nt — 1 



nr — 1 



nt — 1 



ki = ^ Pij, k H = Pn~f, k G = P G7 , k NS = ^ e 2 I P iy — - — ^ 
i=i i=i \ n f j=i 



P 



31 



(2-8) 
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Note that the moments of the splitting functions Pij are related to the anomalous dimensions 
of operators "f n (g) and the coefficient function C % n satisfies the following mass-independent 
RGE: 

'Q 2 m 2 



C3 A^>'j?>W),°)=0, (2-9) 



where r y m (g) is the anomalous dimension for the mass operator. The solution to this equation 
is given by 



ri (Q 2 m " i 2, \ \ T [~ 9ip2) , ln(g,a) 



'j 



m 



where the anomalous dimension "f n (g, a) is a 5 x 5 matrix and is given by 



7n(g,a) = 



In 

K n 0, 



In = 



(ill. 






\ 


llH 


Thh 


Tgh 





ll G 


lie 


1GG 





\ 








Tns) 



where 



(2-10) 



(2-11) 



K n = {KtK n H ,K^K n NS ) . 



(2-12) 



which describes the mixing between hadronic operators and the photon operator. 
2.2. Master formula for the moment 

We can summarize our master formula for the n-th moment of the virtual photon struc- 
ture functions in the case where a heavy quark exists as 



M](n,Q 2 ,P 2 )= 

i,j=f,H,G,NS,-y 



4.(1.^,^ 



Texp 



^ dg ln(g,a) 

g{Q 2 ) M 



^2,n\^—gT- 



(2-13) 



where ip (NS) is the flavor singlet (non-singlet) quark operators for the rif — 1 massless 
quarks. H stands for the heavy quark and fn(Q 2 ) is the running mass evaluated at Q 2 . A\ is 
the operator matrix element renormalized at /i 2 = P 2 , while C\ n is the coefficient function 
renormalized at ji 2 = Q 2 . Since A l n , C\ n are the solutions of the renormalization group 
equation, they depend on the running masses fn{Q 2 ) and m(P 2 ). We can decompose the 
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moments Eq. (12-131) into that for the massless case and the additional term AM^^n, Q 2 , P 2 ) 
due to the mass effects: 



Mj(n,Q 2 ,P 2 ) = Ml(n,Q 2 ,P 2 



massless 



+ zW 2 >,Q 2 ,P 2 



(2-14) 



They are obtained up to NLO by the diagonalization of the anomalous dimension matrix, 
which will be discussed in Appendix |A] 

The master formula for the n-th moment to NLO is given by 



M 2 7 (n,Q 2 ,P 2 ) = / dxx n - 2 F^{x,Q 2 ,P 2 
Jo 



a 1 



4vr 2(3 



47T 



a s (Q 2 



1 



a s (P 2 [ 

a s (Q 2 ^ d " +1 ' 



a s (P 2 



i 

+ c n 



a s (Q 2 



0(a s 



(2-15) 



where a (a s ) is the QED (QCD) coupling constant. The summation index i runs over ±, NS 
corresponding to the three eigenvalues A± and A^r 5 of the one-loop anomalous dimension 
matrix in the massless case. While for the present case with a heavy flavor, i runs over 
if), ±, NS for the four eigenvalues of the one- loop anomalous dimension matrix which 
turns out to be X 7 ^, A±, X% s , where we have = X^ s (See Appendix lAl . For the notation 
of the renormalization group parameters we refer to Ref. [T2l) . The LO coefficients the 



NLO coefficients A^B^C n are, 



f>n _ zf{0) pn/i(0) 



A? 



1 



A" - A? + 2ft 

j j 



-Cf ] — - K {0) P?C { ° } ^ 1 ~ ^ 



'2,71 



, ft-(l) pn^ (0) J_ _ no 2 W pn r (0) 



B? 



pnx/ 1 ) pn 
F i 7n P/ c ( ) 



o 



1 j_ fc(0) pnMV 1 
1 + df n 1 2 ' n l + d? 



K-(0) pnr>(0) & dj 



(2-16) 



where K$ [K^') is the 1-loop (2-loop) photon-parton mixing anomalous dimension, Pj n 's 
are projection operators, Cf^ (C^n) * s the tree-level (1-loop) coefficient function and ft 
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{(3i) is the 1-loop (2-loop) beta function, d% = X™/2(3 . = {a/Ai^A^ is the 1-loop 

operator matrix element and C^n ^ s the coefficient function of the photon operator. (See 
the Ref. IT2"j) for details). We can reorganize the summation over i = NS for the case 
with a heavy flavor, into that for i — ±, NS since the eigenvalue A^; gives rise to the same 
exponents in Eq. f 1 2 - 1 5 [) . as A^ = A^ 5 . The additional terms arising from the variation of 
OME and the coefficient functions due to the heavy quark effects are, 



AM](n,Q 2 ,P 2 ,m 2 ) = [ dxx n ~ 2 AF^ (x , Q 2 , P 

Jo 



2 2\ 

,m ) 



a 1 



L i=±,NS 



+ £ ^ 

i=±,NS 



1 - 



« s (p 2 
«,(g 2 ) 



(2-17) 



where AA™, AB™, AC n are the deviations from the massless case due to the heavy quark 
effects. Note that the heavy quark effects do not change the LO coefficients £f (See Appendix 
IA~|) . This can be derived by an alternative method (See Appendix[B]). The explicit expressions 



of AA?, AB?, AC n are 



A A n 

AAl 

AK 1 

ABl 



AC n 



-e 2 H (-12(3 )AAt G (e 2 H -(e 2 ) nf : 

1 0,n 

-ie 2 ,(-12/3 )(e 2 )„ / ^ G ^ 
n 2 + n + 2 9 . 9 



as 



24 



24 



n(n + l)(n + 2) 

n 2 + n + 2 
n(n + l)(n + 2) 



e 2 H (e 2 ), 



6 )n ^l + r/ 
1 



-AB 



NS 



0,n 



XI -X 



xi 



f l + dl 
AB. 



0,n 

7 G V 



XI - A^ 



G 



-l2{3 e 2 H [AB n G + AAl G 
n f 



(2-18) 



where e# is the heavy quark charge, 



+ |, for SU (2) l up- type quark, 
-i, for SU{2)l down-type quark . 



(2-19) 



Here we have adopted the notation of Ref. [6]) for the anomalous dimensions and coefficient 
functions in the MS scheme^ and AA^ ] = 6((e 2 ),0, (e 4 ) - (e 2 ) 2 )AA^ G . 
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Thus there is no change for the moment at LO level. This is explained by the following 
discussion. In terms of our approach, the heavy quark effect is included in the RGE inputs; 
OME and the coefficient functions. The evolution factor is the same as the massless case and 
there is no physical difference which distinguishes the quarks except for the electric charge in 
the RGE inputs at this order (LO). This results is also justified by the explicit calculation for 
the LO moment. So, the heavy quark effect is occurred at NLO level. Therefore the higher 
order(more than NLO) calculation is essential in the case of including the heavy quark effects 
with massless calculations for the moments. The variation of coefficients AA™, AB™, AC n 
can be obtained by considering the possible variation of the coefficient functions and the 
operator matrix element at NLO. 

Here we confine ourselves to the case of the limit: ^Iqcd ^ P 2 ^ m2 Q 2 ■ In this limit 
we have 



n 2 + n + 2 , m 2 1 1 
In 1 

n(n + l)(n + 2) P 2 n n 2 

4 4 n 2 + n + 2 



n 



(n+1) 2 (n + 2) 2 n(n + l)(n + 2) j 



(2-20) 



AB"l = 0, AB n G = 0, AB™ = 2AB%/n f = 0. (2-21) 

The quark (gluon) coefficient functions B^ (Bq) are obtained by taking the difference be- 
tween photon-parton amplitude and the operator matrix elements.^'^ In the large mass 
limit P 2 <C m 2 , the deviations from the massless case are the same for the photon-parton 
amplitudes and the operator matrix elements. Hence we have AB^ = ABq = 0. We will 
discuss the details elsewhere.^* We also note here for the longitudinal structure function 
F2(x,Q 2 , P 2 ) we do not have heavy quark mass effects to the LO (0(a)) and is given by 
the same formula as the massless case. 

§3. Numerical analysis of F^ ff (x, Q 2 , P 2 ) 

The virtual photon structure functions are recovered from the moments by the inverse 
Mellin transformation. In this section we examine the heavy quark mass effects on the 
effective photon structure function^ F^ s (x,Q 2 , P 2 ) defined as 

F: s (x, Q 2 , P 2 ) = F2(x, Q 2 , P 2 ) + fax, Q 2 , P 2 ) . (3-1) 

We evaluate F^ s up to the NLO and compare our theoretical predictions with the existing 
experimental data from PLUTO Collaboration^ and L3 Collaboration.^ For the PLUTO 
(L3) data, we have Q 2 = 5 (120) GeV 2 and P 2 = 0.35 (3.7) GeV 2 . Therefore, we assume that 
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the active flavors are u, d, s (massless) plus c (heavy) for the case of PLUTO and u, d, s, c 
(massless) plus b (heavy) for L3. 

We plot the experimental data from PLUTO group in Fig. [2] and those from L3 group 
in Fig [31 together with our theoretical predictions. We also show the curves of the NLO 
predictions when active quarks are all massless. We use the QCD running quark mass 
m(P 2 ) which is valid up to the NLO 26 ' and we adopt the following values of the quark 
masses as inputs,^ 1 

m c = 1.3GeV (for PLUTO), (3-2) 
m b = 4.2GeV (for L3). (3-3) 

In general the heavy quark mass has an effect of reducing the photon structure functions 
in magnitude. This feature is explained by the suppression of the heavy quark production 
rate due to the existence of their masses. Heavy quark mass effects appear at larger x 
region. Due to the kinematical constraint for the heavy quark production (p + q) 2 > 4m 2 , 
the contribution of heavy quark to the structure functions exists below x max = — and, 

1+ ~Q T 

therefore, the difference between the massless and the massive cases emerges above x max . 
This kinematical "threshold" effect is not clearly seen in our analysis since we adopted the 
framework based on the OPE and took into account only the leading twist-2 operators. 
But still we see that the difference between the massless and the massive cases becomes 
bigger at large x (see Fig J2] and FigJH below). It is also noted that the heavy quark mass 
effects are sensitive to the electric charge of the relevant quark. Since the photon structure 
functions depend on the quark-charge factors (e 2 ) and (e 4 ), the up-type heavy quark gives 
larger contribution to the photon structure functions than the down-type quark. 

In the case of PLUTO data, there is no justification of our approximation for the limit 
Q 2 ^> m 2 . But we find in Fig. [2] that the predicted curve with mass effects shows the trend 
of reducing the "over-estimated" massless QCD calculation, especially at larger x region, 
and appears to be closer to the experimental data. For L3 data, the hierarchical condition 
P 2 <C m 2 <C Q 2 is satisfied. Although the experimental error bars are rather large, we find 
in Fig. [3] that theoretical curves, both massive and massless cases, are roughly consistent 
with the data, except for the larger x region. For the L3 region, the heavy quark mass effects 
are almost negligible since the b has a charge —1/3. 

Finally we present, in Fig. HJ our prediction for the case Q 2 = 30GeV 2 with P 2 = 
0.35GeV 2 , as an illustration when massive charm quark is relevant. The condition P 2 <C 
m 2 <C Q 2 is satisfied. Although there is no experimental data corresponding to this case, we 
find that the heavy quark mass effects are sizable as can be seen in the figure. 
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Fig. 2. F 7 to NLO in QCD and PLUTO data. n f = 4, Q 2 = 5GeV 2 , P 2 = 0.35GeV 2 , 
0.43. 
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Fig. 3. F c 7 ff to NLO in QCD and L3 data. n f = 5, Q 2 = 120GeV 2 , P 2 = 3.7GeV 2 ,x max = 0.63. 

54. Conclusions 



We have investigated the heavy quark mass effects in the virtual photon structure function 
based on the parton picture as well as on the operator product expansion. We have derived 
the master formula for the additional contributions at NLO to the moments due to mass 
effects. The heavy quark mass effect does not change the LO moments of the photon structure 
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functions but it changes NLO moments. We applied this formalism to the phenomenological 
analysis of F^ s (x, Q 2 , P 2 ). We confronted the theoretical QCD prediction to NLO including 
the heavy quark mass effects with the existing experimental data. 

For the kinematical region of the PLUTO data we assumed the 3 quarks (u,d,s) are 
massless, while the charm quark (c) is the heavy quark. In the region for the L3 data we 
took the 4 quarks (u,d,s,c) to be massless, and the bottom (b) quark to be the heavy quark. 
For the PLUTO region, there exists a sizable heavy quark effect at larger x regime. Although 
our approximation assuming m 2 <C Q 2 is not immediately applicable for the PLUTO region, 
the theoretical prediction shows a right trend of describing the experimental data. For the L3 
region, the heavy quark mass effects are almost negligible since the b has a charge —1/3. It 
is somewhat remarkable that the theoretical predictions for the both cases are consitent with 
experimental data as for the total normalization. We only have one adjustable parameter, 
tIqcd, which we took 0.2 GeV. 

It would be interesting to investigate the charm and bottom quark mass effects at the fu- 
ture SUPER-B experiments.^ If the center of mass energy of the future linear collider (ILC) 
is enough to produce the top quarks, its mass effects would be important for the measure- 
ments of the virtual photon structure functions in view of the charge-factor enhancement. 
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Appendix A 

Diagonalization of the anomalous dimension matrix 



We briefly describe the diagonalization of the anomalous dimension matrix ( 12-111) . Ex- 
panding the anomalous dimension in a power series of the coupling constant g: 



g- 



16tt 



A(0) + 
2 in 1 
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16tt 



2\2 In 1 



and decomposing the one-loop anomalous dimension matrix ^i "* as a block matrix 



7 



(o) 



7f 








0,n 

7tvs 



with the 7n expressed in terms of one-loop anomalous dimensions^ 1 

0,n 

U 





Jy(°) 
I n 



0,n n U.ri \ 

7^ 7 GV , ^ 

0,n 0,n 

7^ 7gi/i 

n/— 1 0,n 1 0,n 0,n i 

n f 7^g ^77^g 7gg / 



we get the eigenvalues of the above 3x3 matrix, X n = X 7 ^, A™ , A™ given by 



Introducing the projection operators we can write down 7„ as 



in 



where the projection operators are given as 



nj — 1 rif — 1 



















A± - A™ 



1 / 0,n 



/ n f~ l ^-°' n X" 

n S- 1 (^,n yi\ J_( ' 

nj— 1 0,n 
~^7,/>G 



a;) 
a;) 



0,n 

7gV 

0,n 

7gV> 



\ 



J_ 0,n 



7gg - A T / 



(Al) 



(A2) 



(A3) 



\ n 0,n \ n J 0,n , 0,n i fi 0~n 0,n \ o | * CLn 0,n 1 / \ a\ 

K = 1^ > A ± = o 1 7gg y 17^ - 7gg) + 4 7^g7gV f ■ ( A4 ) 



(A5) 



(A6) 



(A7) 
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Including flavor non-singlet anomalous dimension into the 4x4 matrix ( 1A-2I) . the 
projection operators are extended to 



In 



E 



(A-8) 



where 



/ 



p" 



nj-X 
n f 





1_ 

«/ 

«/ 




pn 



AJ - Al 





n/-l /. O,: 



o \ 






(A-9) 



o 7 



rv-1/ 0,n 



\ /■i/j'i/j 
n./— 1 _,0,n 



A") 



ri f ^ipG 



0,n 

7& - * 







\ 






(A-10) 



J 
















































1 



(A-ll) 



/ 



Now let us see that the leading order coefficients Hi are not changed in the presence of the 
heavy quark mass effects. With the four eigenvalues, the leading-order coefficients Li turn 
out to be 

n 2 + n + 2 



r " 



CI 



24 



n 



f 



24 



n{n + l)(n + 2) rif 
n 2 + n + 2 1 



(A-12) 



n(n + l)(n + 2) AJ - A™ n/ 



{(nz-lX^-i + ^y 



1 



T 'l + d% 



nri 
>~NS 



24 



n 2 + n + 2 
n(n + l)(n + 2) 



(n f -l)((e% f ^-e 2 H ) 2 



1 4- <7 n 

1 + "TVS 



(A-13) 



where we have denoted the LO coefficients with a hat for the massive case in order to 
distinguish them from those for the massless case. Hence the leading-order QCD result for 
the moment of F% is given by 



dxx n - 2 F2(x,Q 2 ,P 2 



14 



a 1 



47T 



1 - 



« s (p 2 ) 



Here we note that 



Ay, = X N s = 7$J> ^ = rf ivs = 7jj/ 2 A) • 



So we find the sum of and £^ 5 : 

^ + X,JVS " n(n + l)(n + 2)l + cE 



y, 



x 



^((e 2 )n^i - e^) 2 , +(n, - l)(e 4 > B/ -i - (n, - iXe 2 ) 2 ^ 



Now let us remind the following relations: 

nj-l 



then we have 



n 2 + n + 2 



rn , rn _ 

L ^ + L NS-^ n{n+1){n + 2)1 + d 



24 n2 + n + 2 !_JyV-l/Ve 



2 2 



1 



e 4 H 



n f 



i> i=i 



We also find £± = C±. So to the leading-order we get 



f dxx n - 2 F](x,Q 2 ,P 2 ) , 
Jo 



1 (*s(Q 2 )\ 

a s (P 2 )J 



2\\ 



which is nothing but the expression for the n / light-flavor case as we expected, 
the above equation, we have used the fact: 



n f -l 



nt — l 



i=i 



i=l 



i=l 



i=l 
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The above result means that there is no difference between the case with n/ — 1 light-flavor 
plus one heavy-flavor and the one with rif light-flavors to the leading-order in QCD. For 
the NLO coefficients AA™, AB™, AC n we perform the similar analysis and get the additional 
contributions given in Eq. (12- 181) . 

Appendix B 

Master formula for the moments in the parton picture 



There is an alternative method to derive the moment sum rule Eq. (I2-17I) which is based 
on the parton picture.^ Consider the n-th moment of the virtual photon structure function 
F^lx, Q 2 , P 2 ) in the case where rif — 1 light quarks and one heavy quark are present. The 
n-th moments of F% is given by 

M 2 » = ql{n)C L {n) + ql{n)C H {n) + G\n)C G {n) + ql s (n)C NS (n) + C\n) , 

(B-l) 

where we have suppressed, for simplicity, the Q 2 as well as P 2 dependence of the moments of 
the structure function, the parton distributions and the coefficient functions. qlr NS \ denotes 
the flavor singlet (non-singlet) quark parton distribution function for the rif — 1 flavors as 
defined in Eq. (12-41) . and G 7 is the gluon parton distribution function. C % {% = L, H, G, NS, 7) 
are the coefficient functions for parton i-type in the virtual photon, 

C L (n) = (e 2 )^ (l + ^-B;) , (B-2) 



C NS (n) = 1 + £b» (B-3) 



47T 

47T 



C H (n) = C L + e 2 H ^AB;, (B-4) 

C G (n) = (e%^B%, (B-5) 

(P{n) = 2/3 {5^ + 3ejj(B^ + AB^)} . (B-6) 

Putting all the above quantities together and noting that the heavy quark distribution differs 
from the light-flavor distribution by an extra contribution Aq n f(n): 

q H (n) = q n f (n) + Aq n f (n) , (B-7) 

we obtain the moment which includes the heavy quark mass effects as 

W{n) = W[n) | m=0 +e 2 H Aq n f(n) + 6(3 e%AB^ 

r\ ry j 

+e 2 H^q nf (n)AB- + e\^G\n)-ABl . (B-8) 
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Here we note that AB™ = lAB^jn]. Denoting r = a s (Q 2 )/a s (P 2 ) we have 



where 



WW = -^^_ ^(n) + _^(„) , (B.9) 
^s(")/ A" = - r*8) + , (B-10) 

A/2(n)/ A = ^ n ( ! - r "" + ) + ^rC 1 - rd ~) + Ad s . ( B " n ) 



7$ - A" 



AAfo = -2/Wi 2)iVS , AA? = -2/?o^H I AA ^ > 
Z0£ s = 2/? ^4 2)JV5 , = 2A)M 2)1 * > (B-12) 



with 



2/3o4 2)7V5 = 12^ G ((e 4 ) - (e 2 ) 2 ) , 

= l2{3 At G (e 2 ) . (B-13) 



The gluon distribution is given by 



G\n)/£- = ^(1 - + ^£-(1 - r** 1 ) , (B-14) 



where 

'0,n 0,n 



L G -\n_ X ni + d n> ^ ~ )n f n (n + l)(n + 2) ' ^'^ 

From these expressions we can derive the (12-171) with (12-181) . 
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